Gapped spin liquid states in a one-dimensional Hubbard model with 
antiferromagnetic exchange interaction 
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We study the phase diagram of a one-dimensional extended Hubbard model with antiferromagnetic 
exchange interaction analytically and numerically. The bosonization and transfer-matrix renormal- 
ization group methods are used in the corresponding coupling regimes. At half-filling, the system is 
a Mott insulator with a finite spin excitation gap if the on-site Coulomb repulsion is fairly smaller 
than the antiferromagnetic exchange J. This Mott-insulator is characterized by the bond-charge- 
density- wave order or spontaneously dimerization. In the weak-coupling regime where the spin- 
charge separation holds approximately, the critical point separating the gapless and gapped spin 
liquid phases is U c ~ J/2. However, as J increases, the spin-charge couplings become important 
and the critical point U c is significantly suppressed and eventually tends to zero as J — > oo. Away 
from half-filling, the charge gap completely collapses but the spin gap persists. 



I. INTRODUCTION 

The notion of spin liquid state was introduced thirty 
years ago when P. Fazekas and P.W. Anderson first 
postulated that due to the frustrated antiferromagnetic 
coupling the Mott insulating state on triangular lattice 
does not break the spin-rotational symmetry [J|. The 
discoveries of high temperature superconductivity and 
other novel correlated many-body phenomena have stim- 
ulated wide interest in spin liquids both theoretically 
and experimentally^. For correlated electrons in two- 
dimensional Cu-0 plane, the most two important energy 
scales are the kinetic energy and the on-site Coulomb re- 
pulsion. However, as far as the magnetic properties of 
a Mott insulator is concerned, the nearest neighbor an- 
tiferromagnetic exchange do play an important roleQ. 
Recently, F.C. Zhang showed that within the Gutzwillar 
approximation, a quantum phase transition from a Mott 
insulator to a gossamer superconductor(which was first 
proposed by R.B. Laughlin4j) at half-filling can take 
place as the on-site Coulomb repulsion is reduced p| . It 
was further anticipated that away from half-filling, this 
gossamer superconducting state can evolve smoothly into 
the resonant-valence-bond (RVB) spin-liquid phase. 

A basic model for investigating the RVB or gossamer 
states is the so-called t-U-J model @, Li 111 III H It is 
an extended Hubbard model by explicitly including an 
antiferromagnetic exchange interaction. While in two di- 
mensions it remains challenging to accurately solve this 
model, it may be instructive to study the corresponding 
one-dimensional system. 

The one dimensional t-U-J model is described by the 
following Hamiltonian 

H = -t (c\ a c i+1 a + fc.c.j + U n i] -n il 

i<7 i 

+J^S i -S i+1 . (1) 



This Hamiltonian is not merely a toy model moti- 
vated by the corresponding two dimensional system, 
it is also relevant to quasi-one-dimensional correlated 
physics. For example, most of the Bechgaard salts (the 
(TMTSF)2X family of quasi-one-dimensional conduc- 
tors) show close proximity of spin-density- wave(SDW), 
spin-Peierls, ferromagnetic, and superconducting phases 
by varying pressures [9j- To investigate the coexistence of 
triplet superconductivity and ferromagnetism in a class 
of quasi-one-dimensional materials, Japaridze et al stud- 
ied the model for ferromagnetic exchange with easy- 
plane anisotropy in the large-bandwidth limit . As a 
byproduct, they predicted a transition to the dimerized 
ordering phase in the case of weak anisotropy by assum- 
ing spin-charge separation in the weak-coupling regime. 
Comparatively, what was less understood is the generic 
feature of the phase diagram in the U-J plane for the 
isotropic antiferromagnetic modelQ. 

At the first glance, the situation with generic U,J>0 
might be trivial. Because in the atomic limit there is 
no explicit frustration between the pair-wise interactions 
U and J, the system is simply a Mott insulator with 
gapless spin excitations. This is in contrast to another 
well-known one-dimensional model system, i.e., the con- 
ventional extended Hubbard model(CEHM) with the on- 
site({7) and nearest neighbor site(U) Coulomb interac- 
tions. In the CEHM, the pair-wise interactions U and V 
are frustrated explicitly, and in the atomic limit, the sys- 
tem is a charge-density-wave(CDW) insulator for large 
V and a SDW insulator for large [/[Tj. However, if the 
on-site repulsion U is reduced, the virtual hopping pro- 
cesses may affect the spin physics in Mott insulators. For 
example, the recent studies have shown that even within 
the CEHM there exists a new phase in between the CDW 
and SDW phases. The new phase is characterized by the 
bond-charge-density-wave(BCDW) order or spontaneous 
dimerization with gapped charge and spin excitations in 
a narrow window extending from the weak to interme- 
diate coupling regimes [T^ lla IT3. Il5l | . Usually, such an 
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insulating phase is caused by the electron-phonon cou- 
pling or explicit frustrations, but now it is caused by the 
electronic correlations only. With this kind of BCDW 
phase in mind, we are going to examine the phase dia- 
gram of the model (1) in order to clarify the competing 
Mott insulator physics in one dimension if the double 
occupation is allowed. 

In the present work, we study the phase diagram of the 
half-filled one-dimensional t-U-J model by using analyt- 
ical and numerical methods, depending on the ratio J/t. 
It is found that the ground state is always a Mott insula- 
tor at half filling and a metal away from half filling. How- 
ever, the spin excitations have been dramatically changed 
by the exchange term. There is a gapped spin liquid 
phase characterized by the BCDW order when the ex- 
change interaction is fairly stronger than the Coulomb 
repulsion at arbitrary filling. The spin gap stems from 
the interplay between the kinetic energy and the anti- 
ferromagnetic exchange interaction. As spin-charge cou- 
plings increase and the bandwidth decreases, the gapped 
spin liquid phase will be suppressed. All these reveal 
a new scenario for the creation of spin gaps in a single 
itinerant electron chain with translational and spin rota- 
tional symmetries and without explicit frustrations |l6j |. 

II. WEAK-COUPLING THEORY: 
BOSONIZATION, RENORMALIZATION GROUP 
EQUATIONS, AND SEMI-CLASSICAL ANALYSIS 

Let us first consider the phase diagram at half fill- 
ing in the weak-coupling regime U/t,J/t -C 1. In 
this re gim e, the bosonization technique may be used 
reliablvfrii fill Il9l l2(i| . We linearize the spectrum and 
pass it to continuum limit by substituting a~ 1 / 2 Cj (J — > 

(2ira)- 1 / 2 Y Jr =+.- eirkFX+iriPrAx) where <Pr,*( x ) are the 
right/lcft-moving bosonic fields. Introducing charge and 
spin bosonic fields, (f> Ctr = (</V,t + Pr,l)/2) <f> s ,r = 
(VV,T — W,4.)/2, respectively, the Hamiltonian density of 
the bosonized model is then given by H = H c + l~ts + 7~(-cs- 
Here, the charge and spin sectors are described by 

H -= T- E i 9 ^,rf + -^(C^ C , + )(C^ C ,_) 

Z7T ' 2-K z a' L 

r=+,- 

--4^cos2^ (2) 

r=+,— 

+^rhi cos2( t ) s (3) 

with <f) v — 4> v ,-i- + 4>u,- for v = c, s respectively. While, 
the spin-charge coupling part is given by 

9cs 

ftcs = - - 9 9 cos 2<p c cos 2(f> s 

- ^(d x (j) c ^ + )(d x (j) c -)cos2(j) s 
+ | £ ^(^0 S ,+)(9 :E S _)cos20 c 
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Here, in order to compare the t-U-J model with the 
CEHM, we follow the weak-coupling g-ology approach 
and adopt the notation of the Ref.[lJ]. In the lowest 
orders of U and J, the scattering matrix elements are 
9i\\ = -aJ/2,g 1± = a(U - J/2), g 2 \\ = aJ/2,g 2± = 
a(U + J/2), ff 3|| = -aJ/2,g 3± = a{U + 3 J/2), and g 4[l = 
aJ/2,g4± — a(U—3J/2). The renormalized velocities and 
Luttingcr couplings of charge and spin sectors are v c = 
2ta+(g i \\+g i ±-g 1 \\)/2w, v s = 2ta+(g^ — ^ 4 _l — 5i||)/ 27r 
and g p = g 2± + g 2 \\ - 9i\\, g a = 92± - 9z\\ + 9i\\, re- 
spectively. The coupling constants g c and g s denote the 
amplitude of the backward and the Umklapp scattering 
of opposite spins, given by g c — g 3 A_,g s = <7i_l respec- 
tively. g P sjgpa (and g C s/gca) come from the backward 
(and Umklapp) scatterings of the parallel/opposite spins, 
respectively, given by g ps = g p<y = g cs = g ca = - J/2 to 
the lowest order in J. 

The low-energy properties of the one dimensional t-U-J 
model in the weak-coupling regime depend on the scaling 
behavior of these coupling constants during the scaling 
a — » ae dl . Tha is, the coupling constants run in terms of 
the following one-loop renormalization group (RG) equa- 
tions 



Xdg p /dl = 


+2 5c 2 + 


9 2 cs + 9sg ps 


(5) 


Xdg c /dl — 


+^9 P 9c 


9s9cs 9cs9ps 


(6) 


Xdg s /dl = 


-2g 2 s - 


9c9cs - 9c S 


(7) 


Xdg cs /dl = 




- 2g P 9cs - 4g s g cs 






-^9c9s 


^QcQps ^Qcsfjps 


(8) 


Xdg ps /dl = 


-^9 P s H 


- 2g p g s - Ag c g cs 






" 


- ^9s9 P s 


(9) 



where A = Airta. Notice that the SU(2) symmetry in the 
spin sector ensures g a = g s , g cs = g ca , and g ps = g pa . 
From these RG equations, one finds that g p ,g c and g s are 
marginal with the scaling dimension 2, while g cs and g ps 
are irrelevant with the scaling dimension 4. So, we may 
first neglect Tl cs , assuming spin-charge separation in the 
weak-coupling regime. In this case, one obtains the RG 
flows of g p , g c from Xdg p /dl = 2(g c ) 2 , Xdg c /dl = 2g p g c 
and g s from Xdg s /dl — —2{g s ) 2 , respectively. Since both 
U and J are positive and g p (0) = U + 3 J/2 > 0, g c 
is always relevant and flows to the strong-coupling fixed 
point g c (l) — * oo. This indicates the existence of the 
charge excitation gap at half-filling. The charge gap is 
given approximately by A c « t\g c / X\ x ^ 2gp when g c <C A. 
For the spin excitations, g s will flow to if g s > 0, or to 
— oo if g s < [2(j. The latter case occurs when J > 2U 
and in this case a spin gap is expected to open. The spin 
gap is given approximately by A s ss iexp(A/2g s ) when 
< -g s -C A. 

Now, we examine how the spin-charge coupling part 
TL CS affects the phase boundary. As U and J increase, 
the g cs coupling first becomes less irrelevant in terms of 
the one-loop RG equations JSJ-JSJ- However, because 
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g c grows with increasing I and dominates over the other 
couplings, the charge excitations are always gapful. Be- 
low the charge-gap energy scale the 4> c field is locked at 
(cos 20 c ) ~ (A c /t) 2< - 1 ~ 9 "/ x \ Thus, we may re-scale the 
cosine term in spin sector, by introducing g* = g s — g cs . 
The one-loop RG equation for g* is obtained simply by 
\dg*(l)/dl — —2(g*(l)) 2 . This implies that the spin-gap 
transition line is shifted to g* = 0. To obtain the phase 
transition line, we follow the same strategy introduced in 
the Ref.[l4j and solve eqs. ©-©numerically by looking 
at which of the couplings g ci g s and g cs becomes rele- 
vant, In our case, g c grows with increasing I faster than 
the others. We stop the integration once g c reaches 1 
and calculate g* — g s — g C s- The positive (negative) g* 
leads to the spin gapless(gapful) state. Here, the vertex 
corrections to the scattering matrices are not included. 
Although they are indeed crucial for the presence of the 
narrow BCDW phase in the CEHM0, their effects are 
limited in the weak-coupling regime only and may be neg- 
ligible in the t-U-J model. The phase boundary obtained 
by this way is shown by the solid line in Fig.l. 
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FIG. 1: The ground-state phase diagram of the t-U-J model 
at half filling. The solid line is the phase boundary determined 
by numerically solving the RG equations which are valid in 
the weak-coupling regime. From the intermediate- to strong- 
coupling regimes, the phase boundary is continued by the 
dashed line which is qualitatively supported by the numerical 
data of finite-size small cluster. The low energy correlations 
are dominated by the SDW and BCDW orders in the gapless 
and gapped spin liquid phases, respectively. 

To characterize the gapped spin liquid phase, 
we examine via a quasi-classical analysis for four 
kinds of the order parameters, the SDW, CDW, 
BCDW, and bond-SDW(BSDW) parameters, defined 
by Osdw = - Ocdw = 

ObCDW = J2a( C j,(T C j+l,<7 + h - C -)' and £>BSDW = 

{— iy(Cj ^Cj+irf — ct ^Cj+ij + h.c). Upon bosonization, 
they are written in terms of bosonic fields as Ocdw{x) oc 
sin 4> c cos <j) s , S dw{x) oc cos 4> c sin 4> s , Obcdw(x) oc 
cos 4> c cos 4>s , and Obsdw( x ) oc sin </> c sin ^> s respectively. 
Neglecting the spatial variations of the fields we focus on 
the following effective potential 

V e ff{<t>cA») = ~9c cos 2(f) c + g s cos 2(j> s 

—g cs cos2(j) s cos2(j) c . (10) 

The couplings g are all effective ones obtained by in- 
tegrating out high-energy degrees of freedom. As the 



charge gap exists everywhere at half-filling and is al- 
ways large than the spin gap(<jr c increases faster than g s 
), both the BSDW and CDW orders are unfavorable in 
the ground state. So in the insulating phase, the SDW 
and BCDW orders compete each other, with the effec- 
tive potential Vsdw = ~9c - 9s + 9cs and Vbcdw 

q, + g s — g cs respectively. It is the SDW(or BCDW) 
which dominates in the ground state for g cs < g s { or 
g cs > 9s)- Therefore, the SDW-BCDW transition line de- 
termined by quasi-classical analysis is the same as that 
of the spin-gap transition line g* = given above. In 
Fig.(l), the deviation of the phase boundary from the 
straight line U = 2 J with increasing J is due to the g cs 
coupling which in turn enhances SDW order. 



III. INTERMEDIATE AND 
STRONG-COUPLING REGIMES 



We believe that the existence of the spin gap is a 
generic feature of the t-U-J model, not limited in the 
weak-coupling regime only. To further clarify the physics 
behind the spin-gap, it is interesting to consider the 
unconstrained t-J model, namely the t-U-J model with 
[7 = 0, in the strong coupling limit J ^> t. In this case, 
one can first bosonize the Heisenberg exchange term and 
taking the hopping term as a perturbation. In the pure 
Heisenberg model, the spin excitation is critical but the 
charge excitation is completely suppressed. Nevertheless, 
the charge field <j) c should be introduced in the bosoniza- 
tion of the Heisenberg term since charge fluctuations are 
now allowed. For the pure Heisenberg model, <p c is pinned 
to 0. However, by introducing a small but finite i-term, 
<j) c will become finite. This will reduce the charge gap and 
at the same time induce effectively a backward scatter- 
ing term in the spin sector. As a result of this backward 
scattering, a spin gap will open. 

The discussions in the preceding section and the pre- 
vious paragraph showed the existence of the spin gap in 
both the weak and strong coupling limits in the uncon- 
strained t-J model. To examine whether this is true in 
the intermediate coupling regime J ~ t, we have calcu- 
lated numerically the zero-field spin susceptibility using 
the transfer-matrix renormalization group (TMRG) at 
half-filling. The TMRG handles directly an infinite lat- 
tice system. It does not have any finite lattice size effect 
and allows a small excitation gap to be accurately deter- 
mined without invoking the finite size scaling[2lj. This 
is an advantage and the reason for us to use this method 
here. The error comes mainly from the truncation of 
basis states in the TMRG iteration. Fig. 2 shows the 
TMRG result of the uniform spin susceptibility for J = t 
and U = 0. In the TMRG calculation, 100 basis states 
were retained and the maximum truncation error is less 
than 10~ 5 . As can be clearly seen from the inset of the 
figure, the spin susceptibility drops exponentially at low 
temperatures. This exponential decay of the susceptibil- 
ity is a direct consequence of the spin gap. 
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the backward scattering of spin excitations and the sub- 
sequent spin gap phase persist. In the limit U — > and 
J — > 0, the spin gap phase exists when the condition 
g s < 0, i.e. J > 2[/cos~ 1 (-|(5), is satisfied. Thus the crit- 
ical exchange constant increases with increasing doping. 
In contrast to the half-filling case, the gapped spin state is 
now metallic. In this Luther-Emery-type state, the CDW 
and singlet pairing (SP) correlations develop asymptoti- 
cally as Ccdw{t) ~ a\r~ 2 +a2 cos(2kpr)r~ Kc , Csp(r) ~ 
a^r~^' Ka , with at being constants of order of 1. Upon 

doping 5 > 0, K c — \J \+ 9 g p ^2^ f~ < 1- ^ course, the 
CDW correlation is the most dominant one in the Luther- 
Emery phase. 



FIG. 2: Zero field magnetic susceptibility of the one dimen- 
sional t-U-J model with J = t = 1 and U = 0. The inset is 
an enlarged low temperature plot. 

So far, we have established the existence of the gapped 
spin liquid phase in the weak/strong-coupling regimes 
and at a special point in the intermediate regime. Now 
we try to see the generic feature of the phase boundary 
in the U-J plane. As shown in Fig.l, the phase boundary 
is very closed to the line J = 2U in the weak-coupling 
regime, but goes below it as J increases, due to the spin- 
charge couplings. We anticipate that the critical line de- 
velops a maximum in the intermediate regime, and then 
approaches zero gradually as J goes to infinity, see the 
dashed line in Fig.l. 

To support this picture and to estimate the cross- 
ing point, we have also studied the spin-gap transi- 
tion by investigating level-crossings of the excitation 
spectra[22l |2j|. This technique was used to determine 
the phase boundaries of the CEHM by NakamurafT^ 
with high accuracy from the numerical data of finite-size 
clusters. The obtained spin-gap transition line is qual- 
itatively in agreement with Fig.l, the crossing point is 
about (1.3,0.35)113. 



IV. AWAY FROM HALF-FILLING: 
WEAK-COUPLING ANALYSIS 

Finally, we briefly address the subsequences after in- 
troducing holes into the system. In the weak-coupling 
regime and away from half filling, the terms result- 
ing from the Umklapp scattering in 7i c and 7i cs dis- 
appear. The other terms are still present but their 
coupling constants become doping dependent. To the 
leading order approximation, g p = a(U + %nCOs(S<5)), 
g s = a(U — 4cos(|5)) and g ps — — a^cos(^<5), where <5 
is the doping concentration. Due to the absence of the 
Umklapp scattering, the charge gap collapses. However, 



V. SUMMARY 



In summary, we have studied the one-dimensional t-U- 
J model analytically and numerically in different coupling 
regimes. At half-filling, the system shows two distinct in- 
sulating phases: a gapless spin liquid phase dominated by 
SDW correlation and a gapped spin liquid phase domi- 
nated by BCDW correlation. The suggested phase dia- 
gram is plotted in Fig.l. In the weak coupling regime, 
the phase boundary line is determined by solving the RG 
equations. This line is anticipated to develop a maximum 
in the intermediate regime and tends to zero at the strong 
coupling limit. If U is larger than a critical value (which 
is estimated about ~ 0.35t, see Ref . pij ) . the gapped spin 
liquid no longer exists, irrespective of the magnitude of J. 
After doping, the gapped spin liquid phase becomes the 
Luther-Emery phase where the charge gap is completely 
suppressed while the spin-gap persists. 

It is remarkable that in the t-U-J model the spin- 
gap behavior resembles to that of the frustrated t-Ji-J2 
model 25, 26] but with a quite different mechanism and 
the charge gap behavior resembles to that of the Hubbard 
model but with an enhanced magnitude. In the t-U-J 
model, the spin gap results from the interplay between 
the antiferromagnetic exchange and the kinetic energy ( 
the constraint of no double occupancy is released). The 
BCDW insulator driven by this mechanism should have 
a large but finite ratio of the charge gap to the spin gap, 
A c /A s . For the band, Kondo and SDW insulators this 
ratio is unit, 1-1.5 and oo, respectively. For the con- 
ventional spin-Peierls systems, this ratio should be very 
large. So, the present BCDW phase may be relevant to a 
class of quasi-one-dimensional Mott-insulators where the 
charge gaps are only several times larger than the spin 
ones. Particularly, we expect that the mechanism may be 
relevant to the Bechgaard salts such as (TMTSF) 2 PF 6 
where a transition from the spin-Peierls phase into the 
SDW phase by changing pressures was observed p^. 
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